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1. 


INTRODUCTION 


This  report  describes  the  results  of  a study  of  Disturbance 

Accommodating  Controller  (DAC)  design.  The  DAC  method  of  design,  using 

a combination  of  waveform-mode  disturbance  modeling  and  state-variable 

control  techniqvies,  was  developed  by  Dr.  C.  D.  Johnson  of  the  University 

1 2 3 A 

of  Alabama  in  Huntsville  ' ’ * . As  a tool  for  design  of  controllers, 
the  DAC  approach  permits  three  primary  modes  of  disturbance  accommoda- 
tion: (1)  cancellation  (absorption)  of  disturbance  effects,  (2)  mini- 

mization of  disturbance  effects,  or  (3)  constructive  utilization  of  the 
disturbances  as  an  aid  in  accomplishing  the  primary  control  task. 

These  disturbance  accommodations  are  realized  in  addition  to  the  usual 
control  efforts  required  to  satisfy  system  performance  requirements 
without  disturbance. 


A large  number  of  papers  have  been  published  which  describe  the 
general  theory  and  procedures  involved  in  the  use  of  DAC  technique. 
However,  only  a few  practical  numerical  examples  have  been  worked  out 
so  far.  The  purpose  of  the  present  report  is  to  carry  out  detailed 
numerical  designs  of  three  DAC  example  problems.  The  results  presented 
in  this  report  were  obtained  using  the  disturbance  cancellation  mode  of 
DAC  design  on  three  linear,  time- invariant  plants.  The  controllers 
obtained  via  this  design  method  are  termed  "Disturbance  Absorber" 
controllers. 

It  is  not  the  intent  of  this  report  to  thoroughly  cover  all  the 
background  theory  involved  in  the  development  of  DAC  design  procedures. 
This  background  theory  can  best  be  obtained  by  reading  the  original 
papers;  see,  for  instance.  References  1-4. 


wwcmiNQ  p»at  bunc-not  ntma> 


2 . SOME  BACKGROUND 


This  report  considers  controlled  systems  (plants)  which  can  be 
described  by  the  state  equations  of  the  form 

x^Ajc+Jiji+F^w 

1=  C x 

where 

x is  the  plant  state  vector, 

11  is  the  plant  control  input  vector, 

w is  the  vector  of  external  disturbance  acting  on  the  plant, 

jr  is  the  system  output  vector,  and 

A,  IJ,  _F,  £ are  appropriate  size,  known  matrices  which  are  assumed 

constant. 

The  class  of  external  disturbances  w(t)  for  which  the  DAC  theory 
is  intended  is  characterized  by  the  presence  of  "waveform  structure"  in 
w(t),  i.e.,  the  functions  w(t)  can  be  described  by  known  differential 
equations  which  the  w(t)  satisfy  "almost  everywhere".  For  purposes  of 
this  study,  the  disturbances  considered  will  be  assumed  to  be  described 
by  the  following  general  set  of  linear  disturbance  state  equations: 

w = H + L x 

z = Dz  + Mx  + 0 
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where 


z is  the  disturbance  ’’state"  vector, 

a is  a sequence  of  randomly  arriving  vector  impulses,  and 

D,  H,  L^,  M are  known,  time-invariant  matrices. 

It  should  be  noted  that  _L  and  M allow  for  "state  dependency"  of  w. 

In  most  practical  applications,  neither  the  complete  set  of  system 
(plant)  state  variables  nor  the  various  components  w^(t)  of  the  disturb- 
ance are  available  for  direct  on-line  measurement.  Therefore,  the 
DAC's  we  seek  to  design  are  restricted  to  operate  only  on  information 
in  the  available  on-line  measurements  of  the  system  outputs,  set-points 
and  servo-commands  and  any  disturbance  components  which  may  actually  be 
available  for  direct  measurement.  One  of  the  most  important  features 
of  the  DAC  technique  is  that  it  can  effectively  handle  disturbances 
which  cannot  be  directly  measured. 

Since  the  idealized  DAC  control  law  is  a function  of  the  real  time 
system  state,  x and  disturbance  state,  z^,  the  required  on-line  data  for 
practical  implementation  of  a DAC  must  be  generated  via  use  of  state 
constructors  (observers)  operating  on  real  time  system  outputs  and 
control  inputs  _u.  Since  the  external  disturbances  w(t)  are  assumed  to 
have  waveform  structure  and  to  be  modeled  by  known  linear  state  models, 
a state  constructor  can  be  designed  to  generate  estimates  jj  of  the 
instantaneous  disturbance  state  z.  In  addition,  that  same  state  con- 
structor can  be  designed  to  produce  estimates  x.  of  the  instantaneous 
system  state  x. 

Thus,  the  designer  can  utilize  such  a composite-type  state  con- 
structor to  implement  DAC  control  laws  in  the  form 

u - f (x,  z,  t). 
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The  first  example  design  problem  we  will  consider  is  a "servo 
problem"  involving  a first-order  controlled  plant  acted  upon  by  a 
random-like  piecewise  constant  external  disturbance,  as  shown  below. 


The  disturbance  w(t)  is  not  known  a priori  and  cannot  be  directly 
measured.  The  servo  command  y is  also  assumed  to  be  piecewise  constant 
and  unknown  a priori,  but  is  assumed  measurable  in  a real  time  fashion. 
The  objective  of  the  control  design  is  to  create  a controller  which 
will  accomplish  high-fidelity  servo-tracking  of  yc  while,  at  the  same 
time,  consistently  counteracting  the  effects  of  w.  For  this  purpose, 
two  separate  DAC  designs  are  to  be  considered:  Case  (A)  - A design 
utilizing  a full-dimensional  composite  state  reconstructor,  and  Case 
(B)  - A design  utilizing  a reduced-order  composite  state  reconstructor. 


A. 


CASE  A;  DAC  DESIGN  WITH  A FULL-DIMENSIONAL  COMPOSITE  STATE 
RECONSTRUCTOR. 


Piecewise  constant  disturbances  w(t)  can  be  accurately  modeled  by 
setting  (see  Reference  1), 

P = P 

where  p is  the  dimension  of  and  p is  the  dimension  of  w. 

L = D = M = 0 

H = I 


therefore 


4^  = Ku  + Kw  - ax^ 
so 


_x  = A_x+ESu.  + Z^.  = - ax^  + Ku  + Kw 
y = C x = Xj^  . 

Hence 

A = [-a] 

= F = [K] 

C = [1J. 


Following  the  ideas  In  Reference  1,  we  agree  to  split  the  control 
11  into  where  u is  that  part  of  the  control  which  will  be  de- 

signed to  counteract  the  disturbance  and  u^R  is  that  part  which  will 
regulate  the  system. 


We  must  first  check  for  the  existence  of  u as  shown  in  Reference  1. 

— c 

The  control  u will  exist  if,  and  only  if. 


P* 


F ; B £ for  some  matrix  _T. 
However,  in  this  particular  example 

B = F = K is  a scalar. 
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Thus,  exists  (T  = 1)  and  can  be  designed  as  (Reference  1) 


= - T w = - w 


(2) 


For  the  design  of  the  full-dimensional  composite  state  reconstructor, 
we  will  use  the  following  form  proposed  in  Reference  1 


A + F _L  + K C 

F H 

(*)- 

*1 

i + 

B 

M + K2  C 

D 

w 

h 

0 

^R 


(3) 


The  gain  matrices  and  Kj  are  designed  by  examining  the  error  dynamics 
which  are  known  to  obey  (Reference  2) 


A + F L + Kj  C 

F H 

M + K C 

D 

(4) 


For  this 


example,  the  preceeding  expression  reduces  to 


(5) 


Let 


% 

A — 

-a  + k^ 

K 

A — 

k2 

0 

*■* 

— 

It  is  desired  that 

£(t) 

= characteristic  matrix  of  the  e dynamics. 


0 rapidly  which  means  the  roots  of  the  char- 
acteristic polynomial  should  be  "large"  negative  numbers.  The  eigen- 

o, 

values  of  A are  found  as  the  roots,  of  the  characteristic  polynomir 


or 


or 

A2  + (a  - k )A  - k2  K = 0 (6) 

Let  the  desired  roots  of  this  polynomial  be  A^,  Thus,  the  desired 

characteristic  equation  is 

(X  - AX)(A  - A2)  = 0 


or 


” A1  ~ X2  = a " kl 


therefore. 


= a + A^  + 


„ *1  X2 
k2  " K 


Choosing  the  roots  A^,  A ^ properly  will  insure  rapid  settling  out  of 
the  error  transients. 


The  composite  state  reconstructor  we  have  designed  thus  becomes 


a + (a  + A^  + Aj) 
~2 


(a  + A^  + A2) 

*r*2 


•a* 


The  next  step  in  DAC  design  is  to  solve  for  u . The  component  u_  is  to 

K “Tv 

be  designed  so  that  ^(t)  follows,  quickly  and  accurately,  any  admissable 
servo  command  ^(t). 

Following  the  suggestion  in  Reference  3,  the  servo  commands  are 
assumed  to  be  generated  by  the  known  dynamical  process 


ic  = G c 


c = E c + y 


where 


£ = "state"  of  the  servo  commands 

G,  JS  = known,  constant  matrices 

JJ  = sequence  of  randomly  arriving  impulses 


Now,  in  our  particular  example,  is  piecewise  constant.  There- 
fore, just  as  in  the  modeling  of  the  disturbance  process,  we  can  choose 

•\ 


*c  = (1> 


C1  = 0£  + ]J  = VI 


r (id 


A necessary  condition  for  achieving  theoretically  exact  servo  tracking 
is  that  the  "trackability"  condition  (Reference  4) 


G = C 0 


be  satisfied  for  some  matrix  0_.  In  this  case,  since 
G = [1],  C = [1] 


then 


[1]  = [1]  0 " if  and  only  if  " _0  * [1]  . 

Therefore,  one  can  choose  0_  = [Ij  and  theoretically  exact  tracking  is 
possible. 

Now,  the  servo  tracking  error  is 
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- 

or 

£y=Gc-Cx=C0c-Cx=C  [ 0 c - x ] 

Let 

£ * 0 c - x 

— ss — 

then 

£ = C C 

— y ss 

or 

ey  = (1)  [ (1)  ci  ~ xi  1 = cx  - x1  . 

(12) 

It  is  desired  that  -*■  0 rapidly.  Thus,  we  must  have 

C1  - > 0 

or 

e ° 

SS 

It  can  be  shown  (Reference  4)  that  the  differential  equation 
obeys  is 

which  G 

ss 

ic  ' A £,5  - B Uj  + ( 0 E - A 0 + 0 ) c - B z + 9 |j 

(13) 

where 

S = - B _T  + F 

13 

i 

• 

This  simplifies  to 


£ = - ae  - Kiu.  + a c.  + p . (14) 

ss  ss  — R 1 

Now,  it  is  proposed  in  Reference  4 that  be  chosen  as  a linear  func- 
tion of  x and  c^  of  the  form 

% = s L x + s 2 c . (15) 

However,  in  our  design  we  will  instead  choose 

-^R  = —1  — ss  + — 2 — * (16) 

In  this  case 

uD  = s.  e + s_  c. 

R 1 ss  21 

and 

e = - a£  - K(s,  e + s„  c,)  + ac,  + p 
ss  1 ss  2 1'  1 

or 

£ - - (a  + K sx)  ess  + (a  - Ks2)  + p . 

In  order  that  £gs(t)  "*  0 f or  arbitrary  c^,  the  term  (a  - Ks^c^  must 
be  made  to  be  zero.  Therefore,  s^  = a/K.  This  leaves 

e = -(a+Ks,)e  + p 

1 ss 

It  is  desired  that  £gg(t)  0 rapidly  between  successive  incoming  im- 

I pulses  of  p(t).  Therefore,  we  need  (a  + K s^)  to  be  a large  positive 

number  (say  in  the  range  5-10).  Therefore,  we  will  set 
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a + K = n * 


n - a 
K 


So  we  now  have 

UR  - al  £ss  + S2  C1  -(^TT4)  S.  + (t)Cl 

‘ (^ir3)  (ci  ‘ V +(t)  ci 


-(ih'M’i  • 

Previously,  u£  was  found  to  be  (2) 

u = - z 
c 


(17) 


Thus , 


u = uc  + Ur  = - z + xx  + (^)  Cj 


(IDEAL) 


From  the  plant. 


(18) 


so 


A 

u 


(PRACTICAL)  (19) 

Since  y ^ and  y can  be 
directly  measured. 


In  referring  back  to  the  expression  for  the  state  reconstructor,  let  us 
substitute  the  derived  expression  for  u in  place  of  u. 
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x * (-  a + a + A^  + A2)x  + K z - (a  + A^  + X^)y  + Ku 


= (Xx  + X2)x  + K z - (a  + A + X2)y  + K 


+ GO  y'l 


= (A.^  + X2)x  - (a  + X^  + X2  - a + n)y  + nyc 


= (A.^  + X2)x  - (Xx  + X2  + n)y  + ny£ 


while  z remains  the  same. 


X A1  X2  a . A1  X2 

z = -—^X+-K~y 


A diagram  of  the  system  with  the  state  reconstructor  included  is  shown 
in  Figure  1. 


It  is  interesting  to  examine  the  transfer  function  interpretation 
of  the  DAC  servo  control  scheme  we  have  designed.  For  this  purpose,  we 
will  now  proceed  to  find  the  transfer  functions 


G = ZISI  and  G = • 

C yc(s)  d w(s) 


To  find  Gc,  let  disturbance  input,  w = 0.  Thus,  we  have 


KS (S  ~ A2) 

/-*  ..  s „2  ./\  _1_\ 


S -(A  +A  -a)S  -a^+A^S-aS  +a(X1+X2>S+nS  -n (A^A^S+A^S+A^n 


KS (S  - A - X2) 


s3  -(X1+x2-n)s2  + [a-n(W  + x^A^j  **  + ^1^ 


From  superposition  then,  the  overall  closed-loop  transfer  function  for 


our  DAC  design  is 


KS  [s  - (X1  + X2] 


= r_nJy  f KS  |_S  - (Xt  4-  X2)J  ~l 

l_S+nJ  c |_s3-(X1+x2-n)s2+|a-n(X1+X2)+xix^j  s+x^x2nj 


It  was  desired  to  check  the  dynamic  performance  of  this  DAC  design 
by  programming  it  on  an  analog  computer  (EAI-221).  For  a check  case. 


Thus , 


c S + 8 


, _ 10S(S  + 18) 

d s3  + 26S2  + 230S  + 640 


10S(S  + 18) 

(S  + 5. 365) (S  + 10.318  + j3.584)(S  + 10.318  - j3.584) 


and 


y = G y + G ,W 
J crc  d 

(As  a matcer  of  interest,  the  Bode  magnitude  and  phase  plots  for  G^  are 
shown  in  Figures  2 and  3. 

An  analog  computer  flow  diagram  of  the  above  case  is  shown  in 
Figure  A.  Results  of  the  simulation  for  this  design  are  presented  in 
Figures  7a,  7b,  8 and  9«  Also  shown  is  data  associated  with  the  design 
for  Case  B of  this  example  (Case  B will  be  discussed  next). 

B.  CASE  B:  DESIGN  WITH  REDUCED-ORDER  COMPOSITE  STATE  RECONSTRUCTOR 


For  some  practical  applications,  where  reducing  equipment  require- 
ments and  costs  is  a critical  design  factor,  a lower  order  state  re- 
constructor may  be  desirable.  A procedure  has  been  developed  (Reference 
2)  for  the  design  of  a reduced  order  composite  state  reconstructor.  The 
order  of  the  reconstructor  designed  using  this  method  will  be  (n  + p - m) 
where  n is  the  dimension  of  the  system  state  vector  jc,  p is  the  dimen- 
sion of  the  disturbance  state  vector  jz,  and  m is  the  dimension  of  the 
system  output  vector 


From  Reference  2,  the  plant  state  estimate  x.  and  disturbance  state 
estimate  z_  are  given  by  the  "assembly  equations" 


« -[cT  (ccV1- I12  jjy+ll2i 


1 = t2  2 (I  - E y ) 


(21) 


(22) 


and  J_  is  generated  by  the  dynamical  equation  (filter  prescription) 


23 


m 


Figure  7a.  Output  response  due  to  disturbance  input. 
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Figure  9.  Output  response  to  W(t)  for  given  servo  command  . 
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(23) 


i 


|=(i)  + £/v)i  + jjiy  + fiu  . 

The  design  of  the  composite  state  reconstructor  for  this  case  proceeds 
as  follows  (Reference  2). 


and  J_22  are  ch°sen  as  any  appropriate  size  matrices  satisfying 
the  following  conditions 


[ c 


0 


RANK 


[&  ■ 


+ p - m 


For  this  problem 


[ l o ][»  0 

_22_ 


is  satisfied  for  T^  = 0 and  T22  any  number,  say  1.  Thus, 


112  = l°J  * I22  = w 


where  F > 0 is  a symmetric  matrix  satisfying  the  following  equation. 

P=PPT+PP-PHTHP+£  with  £ = £T  > 0 
In  this  example 

*11  = pn(0)  + (0)pii  - Pn  KKPn  + <*n  = Pu  r2  + qu  = o 


Let 


qn  - 1 * Pu  = 1/K  • 


So 


P = [1/K] 


and 


±m-i  m<K>  = -| 


(112  + I £ ) B = 


1 K 

2 K 


= - K/2 


(C  CT)_1  C = 1 


»T  . ,,T 


1 = (112  + I C ) ( A C"  - c"  ) - (V  + I tf)  I + £ 


- j (-a)  - (-K/2)  (-i)  = * - t 


V 2 4 


Substituting 


for  u,  we  obtain 

6 = _ K g + fa  K /a  - nS~|  K ^ n 

2 |_2  4 V 2 ;j  »*r-5Vc 

Then 


x = [l]y  = y 

(25) 

2 = § +|  y . 

(26) 

The  estimation  errors  associated  with  this  reduced-order  reconstructor 
are  given  by  (Reference  2) 

x-x=T12£=0 


Figure  5.  Closed- 


C(S  + B) 


Ks+f) 


Thus, 




<S  + a,<S,-K|(^)s-  -HS 


2 

S + aS 


(a  - n)S  + (§  - £ - f + f)  K + | S 


KS = KS 

w s2  + (a-a  + n+^s+^r  s2  + 


Utilizing  the  same  values  as  for  the  full-dimensional  reconstructor 


w S2  + 13S  + 40 


IQS 

(S  + 5)  (S  + 8) 


Thus, 


(s  + s)yc  + [(S  + 5)  (S  + 8)j  W 


This  design  was  also  programmed  on  an  analog  machine.  The  flow 
diagram  is  shown  in  Figure  6.  Results  of  actual  analog  runs  are  given 


in  Figures  7a,  7b.  8 and  9. 


A 


k 


C.  DISCUSSION  OF  THE  ANALOG  SIMULATION  RESULTS 

Figures  7a  and  7b  show  the  response  of  the  plant  to  various  dis- 
turbance inputs  with  no  servo  command  present.  The  sinusoidal  type 
disturbance  inputs  were  not  completely  absorbed  by  the  DAC  controller 
but  the  more  piecewise-constant  (step  type)  inputs  were.  This  is  as  it 
should  be  since  the  DAC  was  designed  specifically  for  piecewise-constant 
disturbances. 

Figure  8 indicates  the  fidelity  of  the  DAC  in  following  servo  input 
commands,  yc(t),  with  no  disturbance  input.  Figure  9 shows  plant  re- 
sponses with  both  disturbance  and  servo  command  present. 

In  an  attempt  to  see  if  the  DAC's  developed  here  could  be  related 
to  classical  lag,  lead,  lead-lag  compensators,  the  following  transfer 
functions  were  calculated  (using  the  design  developed  in  Part  3. a): 

(1)  u/y,  and  (2)  u/yc.  The  output  of  the  DAC  would  thus  be 

u(S)  = Gj  (S)  y(S)  + G2(S)  yc  (S) 

where 

G^S)  = u(S)/y(S) 

G2(S)  = u(S)/yc(S)  . 

First  find  the  transfer  function  for  u(S)/y(S).  The  block  diagram  for 
this  development  is 


I 


This  can  be  reduced  to  the  following  form. 


The  lower  branch  can  be  combined  as  follows: 


1 - 


xi  + a2  + n s - ^ - x2  - \1  - x2  - n 
s - - a2  s - ai  - a2 


S - 2(A  + A ) - n 


Thus , 


u(S) 


y(s)  - 


(Ax  A2/K)[S  - 2(A1  -»-  A2) 
S (S  - Aj  - A2) 


- nl 


y(s) 
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I 


Substituting  the  numbers  used  in  this  example  gives 


u(S)  = 1.4  y(S)  - 


± 2<i  v(s)  = Ti  4 - ^sjl  2v: I y(s)1 

+ 18)  ' J LU  S(s  + 18) J n ;l 


1.4(S  + 18S)  - 8S  - 192 

S(S  + 18) 


1.4S  + 25. 2S  - 8S 

S(S  + 18) 


1.4S  + 17. 2S  ~ 192 

S(S  + 18) 


Factoring  the  numerator: 


- 17.2  ± /295.84  + 1075.2  - 17.2  ± 37 

2.8  " 2.8 


Sx  = -27.1 


S2  = 9.9 


u(S)  _ (S  + 27 ■ 1) (S  - 9.9) 
y(S)  = S(S  + 18) 
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This  is  a non-minimum  phase  transfer  function  and  does  not  conform 
to  the  models  of  any  of  the  previously  mentioned  classical  compensators. 
The  Bode  magnitude  sketch  of  u(S)/y(S)  acts  like  two  classical  lag  com- 
pensators in  series. 

Figures  10  and  11  show  the  Bode  plots  for  this  transfer  function. 
The  break  frequencies  at  u'  = 18  and  u)  = 27.1  are  too  close  to  give  no- 
ticeable effects  on  these  plots. 

Next,  find  the  transfer  function  for  u(S)/y  (S) . The  block  diagram 
for  this  development  is 


^ > 

1 


The  lower  branch  can  be  reduced  to 


f a ~|  rvvfl 
Ls  ■ xi  - x2J  L s J 


*2  n/K 

s (s  - - \2) 


So, 


u(S) 


i yc<s> 


[Aj  A2  ti/k  "j 

S(S  * \ " VJ 


yc(s) 
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(RAD/SFC) 


Substituting  values  for  this  case  gives 


u(S)  = 


•8vc(S)  - 


64 


S(S  + 18) 


yc(s)  = 


[j8  " S(S  + 18)^ 


yc(s) 


and 


u(S)  a .8(S)(S  + 18)  - 64  _ . 8S2  + 14. 4S  - 64 
yc(S)  ~ S(S  + 18)  ~ S (S  + 18) 


Factoring  the  numerator. 


- 14.4  ± 1/207 . 36  + 204.8  _ - 14.4  ± 20.3 

1.6  1.6 


gives 


Thus , 


- 21.7 


3.7 


u(S)  = (S  + 21.7) (S  - 3.7) 
yc(S)  S(S  + 18) 


! 


(30) 


Again,  the  result  is  a non-minimum  phase  transfer  function. 
Figures  12  and  13  give  the  Bode  plots  for  this  transfer  function. 
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4.  EXAMPLE  DESIGN  PROBLEM  NO.  2 

For  the  next  example  design,  an  output  regulator  problem  is  con- 
sidered for  a second-order  plant  subject  to  a "step-ramp-acceleration" 

2 

input  disturbance  w(t)  = Cq  + C^t  + C2t  . The  objective  here  is  to 
design  a DAC,  using  a reduced  order  state  reconstructor,  which  will 
counteract  w(t)  and  make  the  plant  output  y(t)  rapidly  settle  to  zero 
while  obeying  the  closed-loop  characteristic  equation 

y + <*2y  + oc^y  = 0 

where  (ot^,  a^)  are  given.  A diagram  of  the  plant  and  controller  is 
shown  below. 


To  begin,  we  have 

X1  = x2 
X2  = u + w . 

Thus,  for  the  general  form 

x = Ax  + Bu  + Fw 


we  obtain 


and,  from 


y = £ x. 

we  have 


Also,  L = £ , M r ()  since  the  disturbance  is  not  dependent  on  the 

plant  state. 

Now,  it  can  be  shown  (Reference  1)  that  the  state  x.(t)  of  the 
closed-loop  plant  is  governed  by  the  equation 

x = [A  + BK]x  + [£  l - B R | FH]e  (33) 

where  £_  is  governed  by  Equation  (4). 


For  this  example  then. 


Solving  for  the  characteristic  polynomial: 


Comparing  this  to  the  closed-loop  characteristic  equation 


S + a2S  + 


we  see  that  we  must  have 


k2  = “ a2 


kl  = “ ai 


therefore 


Ur  = < * = £-  - a2J. 


and  since  3c(t)  is  assumed  unmeasurable,  for  practical  implementation  we 


must  have 


4 = ai  “ “2^ 


To  describe  the  disturbance  in  state-variable  form,  we  must  first 
find  a differential  equation  which  is  satisfied  by 


w(t)  = CQ  + Cxt  + C2t‘ 


Taking  the  laplace  transform, 


C_  C.  2C_  CrS  + C.S  + 2C, 

0.1.2  0 1 2 

w(s)  = -g-  + — + — = 3 

b SJ  SJ 


Thus,  the  characteristic  polynomial  associated  with  w(t)  is 


= 0 
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Therefore,  we  will  choose  the  matrices  in  the  model 

w = H z 

_z  = D z_  + 

3 

such  that  = I)  _z  has  the  characteristic  polynomial  X =0  and  H has 
the  general  form 

w = Cq  + C^t  + C^t  • 

For  this  purpose,  let 


and 


or 


Therefore,  let 


Now,  develop  a reduced-dimension  observer  to  generate  x,  z. 


From  input  model,  p = 3. 
From  output  m = 1 . 

From  state  model,  n = 2. 
So  (n  + p - m)  = 4. 


Now, 

«-  |cVtV-!ui  lt  + I12l 


(40) 


1 = i22  (i  - I y) 


(41) 


where 


£ = ( 0 + I f£)i  + i y + £ u (42) 


is  an  n x (n  + p - m)  =2x4  matrix 


T^2  isagx  (n+p-m)  =3x4  matrix 
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where 


Substituting,  we  have 


ai 


rank 


0 0 0 0 

10  0 0 
0 10  0 
0 0 10 
.0  0 0 1 


= n 


Continuing 

— / T T 

T *IT  T +T  T 
-12  \— 12  -12  -22  -22 


T T 
-22  -22 


a - [ c i o ] 


= [ 1 0 0 0 0] 


10  0 0 
0 10  0 
0 0 10 
0 0 0 1 
0 0 0 0 


[1  0 0 0]. 


Now,  from  the  error  dynamics  (Reference  2) 


and 

e_  = (0  + £ H ) £ + T22  o 


To  have  approach  zero  rapidly.  It  is  desired  that  (_£  + £ H ) be 
negative  values.  Using  the  relations  for  V and  H_  we  have 


/ 


V- 


10  0 
0 10 
0 0 1 
0 0 0 


11 

'21 

'31 

'41 


[10  0 0 


"large" 
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— — 

“E  0 0 0" 

\ 

0 10  0 

\ 

0 0 10 
0 0 0 1 
0 0 0 0 

+ 

E 0 0 0 

E 0 0 0 

— _ 

_E,,  0 0 0_ 

Thu  . 


Now,  find  the  roots  of  this  characteristic  matrix 


| (V  + E W)  - X ij  = 


(En-X)l 

0 

0 

E21  -X 

1 

0 

E31  0 

-X 

1 

Z41  0 

0 

1 

>-* 

-X 

1 

0 

(S11  -A) 

0 -x 
0 0 

1 

-X 

-(1) 

SO 


-X)(-X3)  -(1)  (E„,  X2  + E. , + E_,  X) 


A4  - In  A3  - ta  A2  - I31  A - Zn 


We  know  that 


0 - (A  - A1)(X  - A2) (X  - X3)(X  - X4) 


0 = A - (A  + A_  + A-  + A. ) A^  + (A. A_  + A_A.  + A A.  + A, A. 

1 2 3 4 12  34  13  14 


+ A^3  + ^X4^X  - (^^^2^3  + + X2X3X4^X  + 


A1A2A3A4 


Equating  coefficients 


^11  ^X1  + + X3  + 


E21  = X1X2  + X3X4  + X1X3  + X1X4  + X2X3  + X2X4 


E31  - (AxA2A3  + X3A2A4  + AjAjA^  + A2A3A4) 


Z41  m X1X2X3X4 


It  is  desired  that  the  roots  have  "large"  negative  values.  Hence,  let 


A = - 6 


A2  = - 8 


X3  = - 10  + j5 
X4  - - 10  - J5 
In  this  case 

E11  = “ 6 " 8 " 10  ~ 10  = " 34 
-Zn  = 48  + 125  + (60  - 6 j 5)  + (60  + 6j5)  + (80  - 8j5) 

+ (80  + 8j  5)  = 453 

E31  = 48  (-10  + j 5)  + 48  (-  10  - j5)  + (-6) (-  10  + j5)(-10  - j5) 

+ (-8) ( - 10  + j5)(  - 10  - j5) 

- - 480  - 480  - 750  - 1000  = - 2710 

= 48  (125)  = 6000 
41 


-34  1 0 0 

-453  010 

-2710  001 

-6000  000 


+ — 22  — 


(43) 
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Proceeding  with  the  remainder  of  the  problem 


( I12  + I C ) B 


0 1 -34 

00  -453 

00  -2710 

0 0 -6000 


0 0 


0 1 -34  0 

00  -453  0 

00  -2710  0 

0 0 -6000  0 


( c cV1  C = [10] 


§* 


[ 1 0 ] = (1)  [ 1 0 ] = [ 1 0 ] 


T T 

i = ( T , + JL  £ ) ( A c ' lf  - cf  )-(£+!  H)£+i 


.//T  . 

C = Z = 0 


-34  1 

-453  0 

-2710  0 

-6000  0 


0 1 1 
0 0 0 


-34  1 0 0 -34 

-453  010  -453 

-2710  0 0 1 -2710 

-6000  000  -6000 


I 


1 

I 

Having  evaluated  these  integrals,  can  now  be  found  as 

§_  = ( V + £ ) _§.  + X + H 

— - 

-34  1 0 0 

-453  010 

-2710  001 

-6000  000 


-703 

1 

1 + 

-12692 

L + 

0 

-86140 

0 

-204000 

0 

= 0 - 


703 

12692 

86140 

204000 


i = 


-703 

-12692 

-86140 

-204000 


Given  this,  21  can  be  calculated  from 


x = [ cT  ( c cT  )_1  - t12  e ] x + t12  £ 


i=T  (1-1  x) 


'0  1 0 ol 
0 0 10  _§_  - 
_0  0 0 u . 


So  for  the  reconstructor 


X1  = y 


X2  = 34y  + §1 


Zf  = §2  + 453y 


z2  = §3  + 2710y 


z3  = + 6000y 


A diagram  of  this  reduced-order  observer  Is  given  in  Figure  14. 

A.  DISCUSSION  OF  SIMULATION  RESULTS 

This  case  was  simulated  on  a digital  computer.  The  simulation 
program  is  shown  in  Appendix  A along  with  a sample  output.  The  results 
of  a simulation  are  shown  in  Figure  15.  For  the  results  shown,  the 
system  output,  y,  was  subjected  to  an  initial  condition  of  five  units. 
The  disturbance  coefficients  used  are  presented  in  Figure  16  and  the 
total  disturbance  magnitude  in  Figure  17. 


I 
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Figure  14.  Reduced-order  state  reconstructor  design.  Case 


o 

•o 

v£> 


o 

o 

o 


o 


o 


o 

o 
• o 

in 
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Figure  17.  Disturbance  input  magnitude  (W)  for  Case 


It  is  evident  from  Figure  15  that  the  disturbance  input  effects 
are  being  accommodated  nicely. 


5.  EXAMPLE  DESIGN  PROBLEM  NO.  3 

The  last  example  design  problem  is  of  the  set-point  regulator  type 
for  a second-order  plant  with  vector  disturbance  as  shown  below. 


The  set-point  input,  y^p,  is  a given  constant  and  the  disturbance  input 
components,  and  w^»  are  of  the  unknown  piecewise  constant  type.  The 
control  design  objective  here  is  to  make  y ■+  ysp  in  the  face  of  the 
given  disturbance  inputs. 

To  begin,  the  differential  equation  governing  y(t)  must  be  derived. 
From  the  given  plant 


y = \ (46) 


such  that  y(t)  consistently  approaches  ygp  for  all  possible  disturb- 
ances, i.e.,  choose  u to  eliminate  disturbance  forms  and  derivatives 
and  complete  missing  terms  in  the  differential  equation. 

So,  in  this  case,  choose 

u = - wL  - w2  + ugp  - k1  y - k2  y ( 


V = - wx  - w2  + usp  - kx  y - k2  y + wx  + w 


- ki  y - k2  y + uc 


y + k2  y + kx  y = usp  (kx  > 0,  k2  > 0) 


It  is  desired  that  Ugp  be  chosen  such  that 


lim  y (t)  = y 


To  find  a suitable  uqp,  proceed  as  follows:  Let 


GSP  = set-Point  error  = y - y 


70 


then 


Sp  = ~ y 

bSP  = ■ Y = " USP  + k2  y + kl  y 

USP  + k2  bSP)  + kl  (ySP  “ eSP) 

k2  CSP  ” kl  CSP  + kl  ySP  “ USP 
It  is  desired  that  £gp(t)  -*■  0 quickly,  so  let 

USP  = kl  ySP 

then 


ESP  + k2  CSP  + kl  eSP  ° 


Laplace  transforming,  we  obtain 


s + k2  s + kx  = 0 . 


A second  order  system  should  factor  as 
(A  - AX)(X  - X2)  - 0 
thus  giving 


A2  - (A1  + A2)A  + A A = 0 


(53) 


(54) 


(55) 


(56) 


Comparing  Equations  (55)  and  (56),  we  see  that 


X1A2 


I 
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i 


fa 


A 


k2  = ~ + ^2^ 


Hence,  choose  A^  and  A^  to  be  "large"  negative  numbers,  say 


Xi  = -5 


a2  = - 10 


kl  = 50 


k2  = 15 


This  results  in 


USP  * kl  ySP  - 50  ySP 


u = - wL  - w2  + ugp  - kx  y - k2  y 


- wx  - w2  + 50  ygp  - 50  y - 15  y 


= - w.  +50  y - 50  x - 15  x. 


This  u is  the  total  value,  i.e., 


u = u„  + u 
R c 


To  derive  the  full-dimensional  state  reconstructor,  we  start  with 
the  state  equations  for  the  plant  and  disturbances. 
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Plant 


thus 

A'[“°]  ‘ 1’[°]  ; !'["»]  ; £'110’ 
Disturbances 


giving 

° 2 » = i = o ■ 

Now,  proceed  as  in  Case  1 for  the  full-dimensional  state  reconstructor 
(Reference  1). 


The  observer  has  the  form 


For  this  example,  this  reduces  to 


For  the  error  dynamics 


with  the  given  quantities 


Therefore 


e = 


11 


12 


21 


1^22 


10  1 
0 10 
0 0 0 
0 0 0 


(£)  + 


(61) 


To  determine  , jC^,  let 


X = 


11 


12 


21 


.22 


10  1 
0 10 
0 0 0 
0 0 0 


= characteristic  matrix  of  the  e dynamics. 


It  is  desired  that  e(t)  0 rapidly;  therefore,  the  roots  of  the  char- 


acteristic polynomial  should  be  "large"  negative  numbers.  To  find  the 


% 


eigenvalues  of  A,  we  proceed  as  follows: 


lA  - XI  | = 


(kn  -X)  1 0 1 


12 


21 


22 


-X  1 0 

0 -X  0 
0 0 -X 


= 0 


(ku  -xu-x3)  - i (>.2k1?  + k21  X)  - 1 (X2k22)  ^ 0 


X ~ ku  X - k12  X " k22  X ~ k21  X = 0 
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- 


(62) 


A3  - ku.  A2  - (k12  + k22)A  - k21  = 0 
Let  the  roots  be  A.^  A 2,  A3»  Then 

(A  - Al)(A  - A2)(A  - A3)  = 0 

(A  - Al)[A2  - (A2  + A3)A  + A2A3]  = 0 


*3  ~ + A2  + A3)A2  + (A1A2  + A3A3  + A2A3)A  - ^xA2A3  = 0 . (63) 

Comparing  Equations  (62,  63),  we  see  that 


kH  “ " (Xi  + A2  + V * kll  - \ + X3 

- k21  = - A3A2A3  * k21  = A3A2A3 

-ki2  - k22  = a3a2  + a3a3  + a2a3  * k12  + k22  = - (a1a2  + a3a3 

+ a2x3)  . 

Let 


A2  = - 8 + J4 


X3  = - 8 - j4 


Then 


kH  = - 5 - 8 + j4  - 8 - j4  = - 21 


k2]_  = - 5 [SO]  = - 400 


kl2  + k22=  - (40  - 5 j 4 + 40  + 5 j 4 + 80]  = - 160 


k^2  - - 160  - k22 


Let 


22 


80,  then  k^2  = - £50 


Thus 


*2 


[-400 

l -80/ 


Substituting  back  into  the  expression  for  the  reconstructor 


— 

— 

/xl\ 



-21 

1 

0 

1 

-21 

0 

-so 

0 

1 

0 

/ x2 

-so 

y + 

l 

-400 

0 

0 

0 

1 7 

-400 

0 

-so 

0 

0 

0 

\ l! 

-so 

0 

V 

— — 

(64) 


and 


21  + 22  + 2!  y 
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<1 


80  x^  + z^  + 80  y + u 

z^  = - 400  x^  + 400  y 

z2  = - 80  + 80  y 

Thus 

u = 50  ysp  ~ 50  x1  - 15  x2  - 15  z2  - zl  (65) 

is  the  sought  control. 

A block  diagram  of  the  system  with  the  state  reconstructor  is 
shown  in  Figure  18.  This  problem  was  programmed  on  an  EAI-221  analog 
computer  for  checkout  of  the  results.  Table  1 lists  the  programmed 
variables,  the  maximum  values  expected  and  the  scaled  value  for  each. 
Appendix  B shows  the  scaling  process  for  programming  the  analog  and 
Figure  19  gives  the  patching  diagram  for  the  problem. 

A.  DISCUSSION  OF  ANALOG  SIMULATION  RESULTS 

Figure  20  shows  the  plant  response  to  individual  inputs  y , w , 

SP  1 

and  w2«  Response  is  good  to  changes  in  the  setpoint  and  to  changes  in 
w2<  Response  to  w^  is  more  sluggish.  Figure  21  further  emphasizes  the 
sluggish  response  to  but  indicates  that  the  disturbance  is  nulled 
out  when  given  time  between  impulses  in  w^. 

Figure  22  has  output  for  simultaneous  applications  of  and  w2 
with  an  off-zero  setpoint.  Note  chat  when  time  between  impulses  in  the 
disturbances  is  allowed,  the  disturbance  effect  is  nulled  out. 
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Figure 


. 


N >•  N>  N>*  X>  X>  • X>  X>» 


TABLE  1.  SCALED  VARIABLES  FOR  ANALOG  SIMULATION,  PROBLEM  NO.  3 


VARIABLE 

MAXIMUM  VALUE 

SCALED  VARIABLE 

1 

WI 

100 

[ 

1.0  wL  ] 

W2 

100 

[ 

1.0  w2  ] 

ySP 

100 

f 

1*°  ysp  1 

* 

xi 

300 

1 

0.3334  & ] 

A 

X1 

200 

[ 

0.500  xx  ] 

• 

A 

X2 

4000 

[ 

0.025  x 2 } 

A 

X2 

200 

f 

0.500  x2  ] 

X 

Z1 

1500 

[ 

0.06667  zx  ] 

A 

Z1 

120 

f 

0.8334  z1  ] 

X 

z2 

300 

[ 

0.334  z2  ] 

A 

Z2 

100 

[ 

1.0  z2  ] 

y(x2) 

100 

[ 

1.0  y ] 

U 

4000 

[ 

0.025  u ] 

xi 

200 

( 

0.500  x1  ] 

X2 

4000 

f 

0.025  xL  ] 

300 

f 

0.3334  y ] 
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Figure  21.  Output  response  to  individual  w and  w inputs 
with  given  y . * 
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Figure  22.  Output  response  to  simultaneous  applications  of  w and  w 


6. 


CONCLUSIONS 


The  DACs  arrived  at  for  each  example  (using  the  methods  developed 
in  the  references)  performed,  in  each  case,  those  functions  which  they 
were  designed  to  perform.  The  effects  of  unmeasurable  disturbance  in- 
puts were  canceled  out  by  the  portion  of  the  control  which  was  designed 
specifically  to  handle  a given  waveform  mode  disturbance.  When  the  im- 
pulse train  was  not  of  too  high  a frequency,  the  errors  engendered  by 
the  disturbances  were  settled  out  very  well. 
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APPENDIX  A 


This  appendix  contains  a listing  of  the  computer  simulation  of 
Example  No.  2 and  a sample  output.  In  this  program,  the  following 
symbols  apply: 

XI  = output  of  first  integrator  after  u + w 
X2  = Y = system  output 


X5  = §3 

X6  = §. 
A 


PROGRAM  MAIN  74/74  OPT=l 


FTN4. 2+74355  04/21/78  10.59.20 


! 


5 


10 


15 


20 


PAGE  1 

PROGRAM  MAIN ( INPUT,  OUTPUT,  TAPE5=INPUT,  TAPE6=OUTPUT) 

C INTEGRALS  FROM  RUNGK 

COMMON  XI,  X2,  X3,  X4,  X5,  X6 
C DERIVATIVES  TO  RUNGK 

COMMON  XD1 , XD2,  XD3 , XD4 , XD5 , XD6 

COMMON  KUTTA,DT,NX 

J=1 

ALF1=10. 

ALF2=50. 

NX=6 

DT=.01 

XD1=XD2=XD3=XD4=XD5=XD6=0. 

X1=X3=X4=X5=X6=0 . 

X2=5. 

W=U=UC=UR=0. 

WRITE  (6,2000) 

2000  FORMAT  (2X, *TIME*, 12X, *Y*, 15X, *W* , 15X,*U*, /) 

1010  CONTINUE 

IF(TIME.GE.l.)  GO  TO  10 
C0=-1. 

C1=0. 

C2=l. 

GO  TO  50 

10  IF(TIME.GE.2.)  GO  TO  11 
C 0=0. 

Cl=l. 

C2=-l . 

GO  TO  50 

11  IF(TIME.GE. 3.)  GO  TO  12 
C0=1 . 

C1=0. 

C2=0. 

GO  TO  50 

12  IF (TIME .GE. 4. ) GO  TO  13 
C0=0. 

Cl=l. 

C2=-l. 

GO  TO  50 

13  IF(TIME.GE.5.)  GO  TO  14 
C0=-1 . 

Cl-1. 

C2=0 . 

GO  TO  50 

14  C0=1. 

C1=0. 

C2=0. 
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PROGRAM  MAIN  74/74  OPT=l 


FTN4. 2+74355  04/21/78  10.59.20. 


1 


50  CONTINUE 

IF(TIME.GE. 6. ) GO  TO  1000 
DO  100  KUTTA=1 ,4 

50  W=C0+C1*TIME+C2*TIME*TIME 

U=UC+UR 
XD1=W+U 
XD2=X1 

XD3=-703. *X2-34. *X3+U+X4 
55  XD4=-453. *X3-1 2692 . *X2+X5 

XD5=-86140. *X2-2710. *X3+X4 
XD6=-6000. *X3-204000. *X2 
Z1=X4+453.*X2 
UC=-Z1 

60  UR=-ALF1*X2-ALF2* (X3+34 . *X2) 

Y=X2 

GO  TO  (30,60,30,40) ,KUTTA 
30  TIME=TIME+.5*DT 
40  CONTINUE 

65  60  CALL  RUNGK 

100  CONTINUE 

WRITE  (6,2001)TIME,Y,W,U 
2001  FORMAT  (2X,4 (G12 . 6, 3X) 

GO  TO  1010 

70  1000  CONTINUE 

END 


Page  2 
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SUBROUTINE  RUNGK  74/74  OPT=l 


FTN  4.2+74355  04/21/78  10.59.21 


I 


Page  1 

SUBROUTINE  RUNGK 
COMMON  X , DX , KUTTA , DT , NX 
DIMENSION  X(6) ,DX(6) ,XA(6) ,DXA(6) 

GO  TO  (10,30,50, 70) .KUTTA 
5 10  DO  20  1=1, NX 

XA(I)=X(I) 

DXA(I) =DT*DX(I) 

20  X(I)=X(I)+.5*DXA(I) 

RETURN 

10  30  TDT=2. *DT 

HDT=. 5*DT 
DO  40  1=1, NX 
DXA ( I ) =DXA ( I ) +TDT*DX ( I ) 

40  X(I)=XA(I)+HDT*DX(I) 

15  RETURN 

50  DO  60  1=1, NX 
VDT=DT*DX(I) 

DXA ( I )=  DXA ( I ) +2 . * VDT 
60  X(I)=XA(I)+VDT 
20  RETURN 

70  DO  80  1=1, NX 

80  X(I)=XA(I)+(DXA(I)+DT*DX(I))/6. 

RETURN 

END 
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$$ 

$$ 

$$ 
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$$$$$$$ 

$$$$$$$$ 

$$$$$$ 

$$ 

$$ 

$$$$$$$$ 

$$ 

$$ 

TIME  Y W U 


•500000E-01 

.326144 

-.997500 

1736.07 

. 100000E+00 

-2.20212 

-.990000 

1064.53 

.150000 

-2.09702 

-.977500 

278.066 

.200000 

-1.24473 

-.960000 

-31.3571 

.250000 

-.442264 

-.937500 

-109.358 

. 300000 

. 987918E-01 

-.910000 

-101.923 

.350000 

.388633 

-.877500 

-72.7892 

.400000 

.496384 

-.840000 

-44.9627 

.450000 

.490198 

-.797500 

-24.2326 

. 500000 

.421491 

-.750000 

-10.3340 

.550000 

.324981 

-.697500 

-1.55719 

. 600000 

.222706 

-.640000 

3.69061 

.650000 

.127913 

-.577500 

6.55314 

.700000 

.478995E-01 

-.510000 

7.79105 

.750000 

-.  140812E-01 

-.437500 

7.91664 

.800000 

-.575379E-01 

-.360000 

7.29715 

.850000 

-.838194E-01 

-.277500 

6.21361 

.900000 

-.954126E-01 

-.190000 

4.88815 

.950000 

-.953897E-01 

-. 975000E-01 

3.49451 

1.00000 

-.869770E-01 

-. 312639E-12 

2.16203 

1.05000 

-. 732859E-01 

- . 525000E-01 

1.02897 

1.10000 

-. 571629E-01 

-.110000 

.235651 

1.15000 

-. 407288E-01 

-.172500 

-.292992 

1.20000 

-. 254524E-01 

-.240000 

-.610963 

1.25000 

-.122873E-01 

-.312500 

-.750882 

1.30000 

-.175652E-02 

-.390000 

-.742044 

1.35000 

. 597311E-02 

-.472500 

-.615799 

1.40000 

. 1 10129E-01 

-.560000 

-.404780 

1.45000 

. 136710E-01 

-.652500 

-.140298 

1.50000 

.143752E-01 

-.750000 

.150097 

1.55000 

. 1 36042E-01 

-.852500 

.444192 

1.60000 

. 118334E-01 

-.960000 

.725774 

1.65000 

. 949435E-02 

-1.07250 

.984509 

1.70000 

•694935E-02 

-1.19000 

1.21529 

1.75000 

. 447902E-02 

-1.31250 

1.41729 

1.80000 

. 228023E-02 

-1.44000 

1.59283 

1.85000 

•471931E-03 

-1.57250 

1.74638 

1.90000 

-.893733E-03 

-1.71000 

1.88351 
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1.95000 

-.  181761E-02 

2.00000 

-. 233887E-02 

2.05000 

-. 415149E-03 

2.10000 

•435080E-02 

2.15000 

. 718942E-02 

2.20000 

. 725833E-02 

2.25000 

. 543347E-02 

2.30000 

. 280351E-02 

2.35000 

.174617E-03 

2.40000 

-. 199187E-02 

2.45000 

-. 350440E-02 

2.50000 

-.434698E-02 

2.55000 

-.460272E-02 

2. 60000 

-.  440078E-02 

2.65000 

-.388231E-02 

2.70000 

-.  317993E-02 

2.75000 

-.  240644E-02 

2.80000 

-.164973E-02 

2.85000 

-.  971774E-03 

2.90000 

-.410388E-03 

2.95000 

. 173459E-04 

3.00000 

. 310808E-03 

3.05000 

-. 439831E-02 

3.10000 

-.150811E-01 

3.15000 

-. 207983E-01 

3.20000 

-. 198250E-01 

3.25000 

- . 144049E-01 

3.  30000 

-. 721390E-02 

3.35000 

-. 209722E-03 

3.40000 

. 550040E-02 

3.45000 

•947118E-02 

3.50000 

. 116878E-01 

3.55000 

. 123785E-01 

3.60000 

. 118836E-01 

3.65000 

. 105715E-01 

3.70000 

.878648E-02 

3.75000 

. 682119E-02 

3.80000 

. 490318E-02 

3.85000 

. 319274E-02 

3.90000 

. 178744E-02 

3.95000 

. 731072E-03 

4.00000 

. 247644E-04 

4. 05000 

. 100493  E—  0 1 

4.10000 

. 327259E-01 

4.15000 

. 446498E-01 

4.20000 

. 422447E-01 

4.25000 

. 304107E-01 

4.30000 

.149235E-01 

4.35000 

- . 210559E-04 

-1.85250 

2.01015 

-2.00000 

2.13192 

1.00000 

.231780 

1.00000 

-1.88795 

1.00000 

-2.16349 

1.00000 

-1.77639 

1 . 00000 

-1.32481 

1 . 00000 

-.996430 

1.00000 

-.810407 

1.00000 

-.733965 

1.00000 

-.728246 

1 .00000 

-.762274 

1.00000 

-.814672 

1.00000 

-.871784 

1.00000 

-.925403 

1.00000 

-.971005 

1.00000 

-1.00654 

1.00000 

-1.03162 

1.00000 

-1.04694 

1.00000 

-1.05389 

1.00000 

-1.05418 

1.00000 

-1.04964 

-6.25250 

3.64249 

-6.51000 

8.76494 

-6.77250 

9.57097 

-7.04000 

8.85373 

-7.31250 

8.01744 

-7.59000 

7.49761 

-7.87250 

7.33350 

-8.16000 

7.44314 

-8.45250 

7.73122 

-8.75000 

8.12178 

-9.05250 

8.56214 

-9.36000 

9.01839 

-9.67250 

9.47000 

-9.99000 

9.90557 

-10.3125 

10.3199 

-10.6400 

10.7121 

-10.9725 

11.0840 

-11.3100 

11.4391 

-11.6525 

11.7820 

-12.0000 

12.1174 

3.05000 

2.45679 

3.10000 

-7.99785 

3.15000 

-9.16160 

3.20000 

-7.06363 

3.25000 

-4.72115 

3.30000 

-3.06286 

3.35000 

-2.16605 
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A. A0000 

- . 1207A5E-01 

3.A0000 

-1.8A939 

A.A5000 

-. 203173E-01 

3.A5000 

-. 190A85 

A. 50000 

- . 2A7537E-01 

3.50000 

-2.167A5 

A. 55000 

-. 2590A3E-01 

3.55000 

-2.52328 

A. 60000 

- . 2A5235E-01 

3.60000 

-2.89938 

A. 65000 

-. 21A17AE-01 

3.65000 

-3.25192 

A. 70000 

-. 1733A3E-01 

3.70000 

-3.55690 

A.  75000 

-. 129050E-01 

3.75000 

-3.80376 

A. 80000 

-.  861629E-02 

3.80000 

-3.9910A 

A. 85000 

A807A2E-02 

3.85000 

-A. 12323 

A. 90000 

-. 168117E-02 

3.90000 

-A.208A7 

A. 95000 

. 675875E-03 

3.95000 

-A.2567A 

5.00000 

. 226859E-02 

A. 00000 

-A.278A5 

5.05000 

-. 27A709E-03 

-1.00000 

-.98001A 

5.10000 

-.  7A08A0E-02 

-1.00000 

2.51706 

5.15000 

-. 112A98E-01 

-1.00000 

2.89988 

5.20000 

-. 105737E-01 

-1.00000 

2.20107 

5.25000 

-.696969E-02 

-1.00000 

1.A28A2 

5.30000 

- .2301 32E-02 

-1.00000 

•8870A7 

5.35000 

. 210199E-02 

-1.00000 

.597906 

5. A0000 

. 55226AE-02 

-1.00000 

. A97356 

5.A5000 

. 770752E-02 

-1.00000 

.51A587 

5.50000 

. 869596E-02 

-1.00000 

. 59AA8A 

5.55000 

. 868376E-02 

-1.00000 

.699791 

5.60000 

. 793059E-02 

-1.00000 

. 8073A9 

5.65000 

. 670221E-02 

-1.00000 

. 9038A6 

5.70000 

. 52372AE-02 

-1.00000 

. 98253A 

5.75000 

. 373028E-02 

-1 . 00000 

1.0A099 

5 . 80000 

. 232582E-02 

-1.00000 

1.07958 

AD-A060  628  ARMY  MISSILE  RESEARCH  AND  DEVELOPMENT  COMMAND  REDSTO— ETC  F/G  9/3 
INVESTIGATION  OF  DISTURBANCE  ACCOMMODATING  CONTROLLER  DESIGN. (U) 

JUL  78  M L MCCOMAN 

UNCLASSIFIED  ORDMI-T-78-65  NL 


END 

DATE 

FUMED 

I 79 

DDC 


NATIONAL  BUREAU  OF  STANOAROS 

MCMOCOPT  RCSOLUTIO N TEST  CHART 


APPENDIX  B 


In  order  to  program  Example  No.  3 on  the  analog  without  having  to 
resort  to  very  high  gains  on  the  amplifiers,  a magnitude  scaling  process 
was  used.  Equations  (47),  (64)  and  (65)  were  programmed  thusly  with  max- 
imum expected  values  being  normalized  to  100  volts.  The  resulting  ex- 
pressions give  the  pot  settings  and  amplifier  gains  for  patching  each  of 
the  equations  on  the  analog. 


x^  = - 21  x^  + $2  + ^2  + 21  y 


[ 0.3334  ^ ] = (0.3334)  ()  21)  [(Q  5)  ^ ^ + (P-.?334)  ^ ] 


+ [l.o  s2  ] + ig-3.3W.C2-1)  tl.o  y] 


[ 0.3334  Xx  ] = - (.14)  100  [ 0.5  5^  ] + (.6668) [0.5  x2  ] 


+ (0.3334)  [1.0  z2  ] + (.7)  10  [ 1.0  y] 


- 80  xx  + Zj  + 80  y + u 


[ 0.025  S2  ] 


(0r025)(-  80)  f , , (0.P25)  f Q ^ g i 

0.5  1 l1  0.8334  1 *JJ  1 J 


+ v [1.0  y ] + [ 0.025  u] 


[ 0.025  i2  ] = - (.4)  10  [ 0.5  ] + (0.3)  [ 0.8334  \ ] 


z2  » - 80  xL  + 80  y 


[ 0.333*  i21  ■ (0.333*),C-_80I  [0.5  ^ , . (0.333*j<801 , 1.0  y ] 


[ 0.3334  \2  ] = - (.5334)  100  [ 0.5  xx  ] + (.2667)  100  [1.0  y] 
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